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We study hysteresis in anti-ferromagnetic random-field Ising model at zero temperature. The 
external field is cycled adiabatically between -oo and oo. Two different distributions of the random- 
field are considered, (i) a uniform distribution of width 2A centered at the origin, and (ii) a Gaussian 
distribution with average value zero and standard deviation cr. In each case the hysteresis loop is 
determined exactly in one dimension and compared with numerical simulations of the model. 



I. INTRODUCTION 

Hysteresis is a non-equilibrium effect commonly observed in systems subjected to a cyclic force It means that 
the response to a changing force depends on the history of the force. In particular, the response in increasing force is 
different from that in decreasing force. This is caused by the delay in responding to the force. Theoretically hysteresis 
should disappear if the force changes sufficiently slowly but this often corresponds to unrealistically long time periods 
such as the life span of an experimentalist. Several complex and disordered systems like permanent magnets show 
hysteresis over the longest practical time scales. Experience with spinglasses and other systems containing quenched 
disorder Q has revealed that the free energy landscape of such systems comprises a large number of local minima 
(metastable states). The number of local minima is thermodynamically large. The barriers between the local minima 
are also large compared with the thermal energy of the system. Consequently, in the absence of a driving field the 
system gets trapped in one of the local minima and is unable to explore the entire phase space over practical time 
scales. In this situation the thermal relaxation time of the system r is much larger than the relaxation time of its 
constituent units (individual spin- flips) as well as the period 27r/w of the cyclic driving field. A useful approximation is 
to assume t to be infinite or equivalently the system to be at absolute zero temperature. This makes the dynamics of 
the system deterministic and more amenable to analytic solutions and simulations without compromising the essential 
physics of the problem. We take the limit T — > before the limit cj — > to obtain nonvanishing hysteresis in the limit 
w = 0. 

In an extensive and pioneering work Sethna et al used the random-field Ising model Q along with the Glauber 
dynamics 0] at zero temperature to study hysteresis in ferromagnets with quenched disorder. They analyzed their 
model using numerical simulations, mean field theory, Wilson's renormalization group and compared it with 
experiments. Their model reproduces several experimentally observed features. These include familiar shapes of 
hysteresis loops, Barkhausen noise and return point memory. Interestingly the model predicts the existence of 
a non-equilibrium critical point on each half of the hysteresis loop. This is based on a Gaussian distribution of the 
random-field with mean value zero and standard deviation cr that plays the role of a tuning parameter in the model. 
The model may be solved exactly in one dimension and on Bethe lattices of a general coordination number z 

EMI- 

Above a lower critical coordination number z [T^.[T3j there is a critical value cFc such that for a < cFc, each half of the 
hysteresis loop has a first order jump in the magnetization at some applied field h. The size of the jump goes to zero 
as cr —>■ CTc from below. If he is the critical field at which the jump vanishes, {/ic, (7c\ is a non-equilibrium critical point 
showing scaling of thermodynamic functions and universality of critical exponents in its vicinity. This is reminiscent 
of equilibrium critical phenomena and a ppe ars to have a fair amount of experimental support in the field of hysteresis 
as well. A generalization of the model |14l4l7l | to n-component [n > 1) classical spins shows the existence of critical 
points in the generalized model as well. The critical exponents of the generalized model are in the universality class 
of the random-field Ising model (n = 1) if the critical point occurs at a non-zero value of magnetization or the applied 
field. This is understandable because a non-zero value of magnetization or the applied field picks a unique direction 
in the system that effectively reduces its symmetry to that of an Ising model. This lends further support to the broad 
agreement between experiments and predictions of RFIM for hysteresis. 

Hysteresis in the anti-ferromagnetic random field Ising model [l^-[2lj has received relatively little attention as 
compared to its ferromagnetic counter part Ej, H^-HSl ■ This is partly due to the difficulty of obtaining 
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analytic solutions in the anti-ferromagnetic case. For the ferromagnetic case, exact expressions for the major and 
minor hysteresis loops have been obtained in one dimension as well as on a Bethe lattice of coordination number z 
for a bounded as well as a Gaussian distribution of the quenched field [l^, HO, [HI • The distribution of Barkhausen 
jumps (avalanches) has also been obtained [l^l- Several other aspects of the ferromagnetic model have been studied 
in the mean field theory as well as on periodic lattices p3 - l27| . In the anti- ferromagnetic case an expression for the 
major loop has been obtained in one dimension in case the quenched field has a uniform distribution of width 
A centered at the origin and A < | J| where J is the anti-ferromagnetic exchange interaction. The purpose of the 
present paper is to extend this result to A > | J| as well. The results presented here are also applicable to unbounded 
distributions of the quenched random-field such as the Gaussian distribution. 

It may appear rather surprising at first sight that there should be any difficulty in solving a one-dimensional Ising 
model at zero temperature. The difficulty arises primarily from the presence of quenched random fields. Problems 
with quenched disorder are difficult to analyze analytically. Besides this the spin-flip dynamics with anti-ferromagnetic 
interactions is more complicated than its ferromagnetic partner. Consider two spin systems of equal size and having 
the same realization of quenched field distribution. Let one system have ferromagnetic nearest neighbor interaction 
J and the other an anti- ferromagnetic interaction — J. In equilibrium, the ground states of the two systems on a bi- 
partite lattice are related to each other by symmetry. Evidently no such relation is available between non-equilibrium 
metastable states of the two systems. An applied field ha increasing adiabatically from ha = —oo to ha = oo takes 
both systems from a stable state with all spins pointing down (i.e. aligned along ha = — oo) to all spins pointing up. 
Although the end points of the trajectory are the same for both systems but the magnetization paths are different. 
In particular the number of metastable states along the two paths arc different. In the ferromagnetic case, spins tend 
to flip up in avalanches and do not flip down in increasing fleld. The anti-ferromagnetic dynamics is marked by the 
absence of avalanches. This is because a spin flipping up at an applied field ha prevents its neighbors from flipping up 
at the same fleld. However, a spin flipping up at ha occasionally causes its neighbor to flip down at ha- This is a kind 
of a reverse avalanche in anti-ferromagnetic dynamics that involves only two spins including the spin that triggers the 
avalanche. The forward avalanches in the ferromagnetic case, and the reverse avalanches in the anti-ferromagnetic case 
provide a mechanism for irreversibility in the two models respectively and give rise to hysteresis. Due to the smaller 
size of reverse avalanches the area of anti-fcrromagnetic hysteresis loop is much smaller than the area of ferromagnetic 
loop. Also the Barkhausen noise on the ferromagnetic hysteresis loop that is caused by large sporadic avalanches is 
nearly absent in the anti-ferromagnetic case. 

The relative difficulty of analyzing anti-ferromagnetic dynamics comes from the fact that it is non-Abelian while 
the ferromagnetic dynamics is Abelian. This means as follows. Consider an unstable system at an applied fleld ha 
such that one or more spins are not aligned along the net fleld at their site. We relax the system till it is stable. 
Relaxing the system means checking each spin and flipping it if it is not aligned along the net fleld at its site. It is 
an iterative process because flipping a spin may reverse the sign of the net fleld at its nearest neighbors. We have to 
continue the relaxation process till each spin in the system is stable. A dynamics is called Abelian if the end result 
of the relaxation process does not depend on the order in which the spins are relaxed. If the result does depend on 
the order in which the spins are relaxed it is called non-Abelian. Consider two nearest neighbor spins which are both 
down but the net field at their sites is positive. If the interaction between the spins is ferromagnetic, the spins can 
be relaxed in any order and the end result would be that both spins are turned up. This is because turning a spin 
up makes the net field at its neighbor even more positive so that the neighbor also has to be turned up. This is not 
the case with anti-ferromagnetic interactions. Turning a spin up decreases the net field at its neighbor and it may 
decrease it below zero so that the neighbor no longer needs to be turned up when relaxed. Thus the end result may 
be one spin up and one down. Which one is up depends on which one was turned up first. The anti-ferromagnetic 
dynamics is therefore non-Abelian. As the stable state at the end of the relaxation process depends on the order 
in which the unstable spins are relaxed, we have to choose a protocol for the order in which the unstable spins are 
relaxed. At every step, we choose to relax the most unstable spin in the system i.e. the one whose flipping would 
lower the energy of the system the most. Locating the most unstable spin at every step of the dynamics is what makes 
the anti-ferromagnetic model more tedious to analyze theoretically as well as numerically. 

Although the aim of the present study is to flnd an analytic solution of a non-equilibrium problem with quenched 
disorder, we may mention some connection with experiments. Relaxation dynamics of any complex statistical system 
belongs to one of two broad categories: (i) where relaxation takes place by avalanches, and (ii) where it proceeds by 
single localized events. The ferromagnetic random-fleld Ising model belongs to the category of avalanches. It explains 
experimental effects such as the Barkhausen noise and the possibility of non-equilibrium critical points. The anti- 
ferromagnetic random-field Ising models belongs to the second category characterized by the absence of avalanches. 
Due to the absence of avalanches, we do not expect small changes in the applied field to cause large changes spanning 
across the system. In other words, we do not expect the response of the system to be critical at any value of the 
applied fleld. This rules out the existence of non-equilibrium critical points in anti-fcrromagncts. Our calculation 
shows that the hysteresis loop of an anti-ferromagnet with relatively small quenched disorder (to be defined in the 
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following) has a wasp-waisted shape i.e. constricted in the middle. In the limit of very small disorder the wasp-waisted 
shape gradually transforms into two hysteresis loops joined by a long and narrow region of almost no hysteresis. For 
much larger disorder the familiar pot-belly shape of ferromagnetic loops is recovered. Thus the anti-ferromagnets 
can exhibit a wide variety of shapes of hysteresis loops and this feature of our model is in general conformity with 
experiments [2^-[13- Anti- ferromagnetic hysteresis loops comprising three loops are also observed in experiments 
[33I [3^ . This too is understandable if our one-dimensional model is extended to lattices with higher coordination 
number. The anti-ferromagnetic model may also apply to other systems that exhibit glassy dynamics [lol . l35l - [38| 
characterized by a single localized events. 

The outline of this paper is as follows. In Section II we define the model for a general distribution of the quenched 
random-field centered at the origin. We focus on two specific cases: a uniform bounded distribution of width 2 A 
and a Gaussian distribution of standard deviation a. If | J| is the magnitude of nearest neighbor anti- ferromagnetic 
interaction, the case A < | J| is easier to treat analytically for reasons to be made clear in the following. The case 
A > I J| is more difficult but comparable to the case of unbounded Gaussian distribution. In Section III, we show 
numerical results for the hysteresis loop in three representative cases: (i) A ~ 0.5| J|, (ii) A = 1.25|J|, and a = .5| J|. 
Case (i) lies in the ambit of an exact solution obtained earlier [13, for A < | j|. In Section IV, we briefly review 
the earlier result because it is needed to proceed to distributions with A > |J| and the Gaussian distribution. The 
formalism and presentation of section IV is slightly different from the earlier version [2^ [2l| on which it is based in 
order to make a smoother transition to the following sections. Sections V, VI and VII treat a uniform distribution of 
arbitrary A and a Gaussian distribution of arbitrary a . These sections contain the main results of this paper. Exact 
expressions for the hysteresis loops arc obtained. These have been superimposed on the simulation results shown in 
figures (l)-(6). The fit between the theory and simulations is quite good as may be expected from an exact solution. 
Indeed the two are indistinguishable on the scale of the figures. The fact that simulations over a relatively small size 
of the system agree with the exact result is due to the super-exponential decay of correlations in this system [sll . 
The agreement between simulation and theory also justifies (albeit post facto) the implicit assumption in our analysis 
that the system is self-averaging. 



We consider non-equilibrium anti-ferromagnetic random-field Ising model in one dimension at zero temperature. 
At each site i (i = 1, 2, 3, ... , N) of a linear lattice, there is an Ising spin = ±1 which interacts with its nearest 
neighbors through an anti-ferromagnetic interaction J (J < 0). A quenched random-field hi as well as a uniform 
externally applied field ha acts on si. The Hamiltonian of the system is. 



II. THE MODEL 




(1) 



We consider two distributions ^(ft.,) of the random-field \hi\: 



(a) A uniform bounded distribution of width 2A centered at the origin. 



ci>(hi) ifhA</i, <A] 



~ 0, otherwise. 



(2) 



(b) A Gaussian distribution with average zero and standard deviation ct. 



1 




(3) 



It is convenient to rewrite H in terms of the net field /, acting on spin s. 




(4) 



The spins {si{t) ~ ±1} obey discrete-time single-spin-flip Glauber dynamics at zero temperature i.e. Si{t + 1) = 
sign fi(t). This means a spin flips only if it lowers its energy. It also assumes that if a spin-flip is allowed, it occurs 
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at a rate F which is much larger than the rate at which the magnetic field ha is varied. Thus all flippable spins relax 
instantly and Si{t + 1) has the same sign as the net local field fi{t) at its site. 

s,(t + 1) = sign f,{t) = sign [J{s,^i{t) + + h, + ha{t)] (5) 

Iterative application of the above dynamics leads to a fixed point state of the system such that Si{t + 1) = Si{t) 
for each spin Si (t) in the system. The condition of adiabatic variation of the applied field or equivalently the instant 
relaxation of spins mentioned previously is implemented by holding the applied field constant until a fixed point is 
reached. The fixed point is a local minimum of the energy (metastable state) of the system. We denote the fixed 
point value of Si{t) by s* and characterize the fixed point state by the magnetization m{ha) per spin, 

■m{ha) ^ j^^s* (6) 

i 

Our aim is to find the magnetization m{ha) of each metastable state visited by the system as the applied field is 
cycled adiabatically from ha = —oo to ha = oo and back to ha = —oo. We start with ha = —oo when each spin is 
necessarily aligned along the applied field i.e. we have a fixed point with m ~ — I. Now we increase ha slowly till some 
spin becomes unstable and needs to be flipped. We flip this spin and check its neighborhood if any more spins have 
to be flipped. If any neighbors are flipped, we check their neighbors if they need to be flipped as well. This process 
is continued till no more spin needs to be flipped i.e. we reach a new fixed point. The applied field is held constant 
during the passage from the old fixed point to the new. This procedure is continued up to ha = +oo when m = 1. 
The magnetization m^(h a) on the upper half of the hysteresis loop when ha is slowly decreased from ha — to 
ha = — oo can be obtained from m{ha) by a symmetry relation my'{ha) = -~'m{—ha). Therefore the calculation of the 
lower half of the hysteresis loop suffices to determine the entire hysteresis loop. 

III. SIMULATIONS 

Computer simulations of the preceding model play a useful role in guiding its analysis and checking the analytic 
results. Normally we used 10"^ spins on a linear lattice with periodic boundary conditions and used 10^ independent 
realizations of the random field distribution to generate the data. The data was binned in 10^ bins in the applicable 
range of the applied field and averaged over different realizations of the field distribution. The simulation results 
shown in figures (l)-(6) and some of the other figures in section VIII were obtained in this way. It took approximately 
four hours on our 3 GHz desktop to generate the data for each figure. When the estimated probability of an event 
was very small, say of the order of 10~^, we performed the corresponding simulation on a larger system, say 10^ spins, 
and a smaller number of independent runs, say 10^. This was to optimize the accuracy of the data and the time taken 
to generate it. In some cases it was more appropriate to have a larger number of bins. In these cases we worked with 
the data obtained from a large system of 10^ spins without binning it. We set J = —1, and as mentioned previously 
we performed simulations for three cases: (i) A = 0.5, (ii) A = 1.25, and (iii) a = 0.5. These values were chosen 
arbitrarily but represent three broad classes of the analytic results. 

A. A = 0.5 

Figure (1) shows the hysteresis loop for A = 0.5. We see that m{ha) = —1 if ha < — 2| J| — A, and m{ha) = 1 if 
ha > 2| J| + A. The magnetization m(ha) rises from -1 to +1 in three steps. We call these steps ramp-I {ha = — 2| J| — A 
to ha = -2| J| + A); ramp-II {ha = -A to ha = +A); and ramp-III {ha = 2\J\ - A to ha = 2\J\ + A). The ramps 
are connected to each other by two plateaus; plateau-I ( ha = — 2| J| + A to ha = —A); and plateau-II {ha = +A 
to ha = 2\J\ — A). On the plateaus, the magnetization remains constant even though the applied field continues to 
increase. Plateaus occur for A < | J| (small disorder), and simulations suggest that magnetization on the plateaus 
is independent of A. Numerically, the magnetization on the plateaus is approximately = —.135 on plateau-I, 
and m^^ = .109 on plateau-II. The qualitative shape of m{ha) is easy to understand. Due to the anti- ferromagnetic 
interaction between nearest neighbors, spins with both neighbors down are the first to turn up in an increasing applied 
field. Such spins turn up on ramp-I. Next are the spins with one neighbor up and one down which turn up on ramp-II. 
Spins with both neighbors up require the largest applied field to turn up, and these turn up on ramp-III. For A < | J|, 
the three ramps are well separated from each other. In other words, no spin with n up neighbors (n = 1,2) can turn 
up in increasing ha until all spins with rt — 1 up neighbors have turned up. On each ramp, the sequence in which 
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the spins turn up is determined by the distribution of the quenehed random field. Spins with large positive quenched 
field turn up before spins with a lower quenched field. The quenched field lies in the range —A to +A. Thus each 
ramp has a width 2A along the axis of the applied field. When a spin turns up on ramp-I. its nearest neighbors 
are placed in a category so that they cannot turn up before ramp-II. Similarly when a spin turns up on ramp-II, 
its nearest neighbor which is down cannot turn up before ramp-III. This is essentially the reason for the absence of 
avalanches in the anti-ferromagnetic RFIM. Occasionally on ramp-II and ramp-III, a spin turning up can cause its 
nearest neighbor which is already up to turn down. We may call this a reverse avalanche of size unity. There are no 
long avalanches as in the ferromagnetic model. If there were no reverse flips at all there would be no hysteresis in 
the model. The smallness of the reverse flips is the reason behind the smallness of the area of the hysteresis loop. In 
order to highlight the separation of the upper and lower halves of the hysteresis loop we have plotted in figure (2) 
the relative separation of the two halves relative to their average value. As the majority of the spins turn up one at 
a time, the calculation of m{ha) becomes essentially a matter of sorting quenched random fields in decreasing order 
on each ramp. The difficulty arises from the fact that a posteriori distribution of random fields on unfiipped spins 
that are next to a fiippcd spin is significantly different from the initial uniform distribution. The main problem is to 
calculate this a posteriori probability distribution of random fields on unfiipped spin sites. 

B. A = 1.25 

Figure (3) shows the hysteresis loop for A = 1.25. We sec that the three ramps comprising the hysteresis loop 
in figure (1) have lost their individual identity. If A > a spin with one neighbor up can turn up on the lower 
hysteresis loop before all spins with both neighbors have turned up. This makes the analysis of the dynamics more 
complex. We shall take it up in Section V and Section VI. Figure (4) is a magnified version of figure (3). It shows 
the relative separation of the two halves from their average value at a given applied field. 

C. cr = 0.5 

Figure (5) shows the hysteresis loop for a = 0.5 for a Gaussian distribution of the random field. The Gaussian 
distribution is an unbounded distribution. Therefore we may not expect the hysteresis loop to comprise of sharp 
ramps and plateaus as in figure (1). However, notice the qualitative similarity between figure (5) for a = 0.5 and 
figure (1) for A = 0.5 sans the sharp edges in figure (1). The general analysis presented in sections V, VI, VII applies 
equally to a uniform distribution with any value of A and a Gaussian distribution for any value of a. Figure (6) is 
a magnified version of figure (5) showing the relative separation of the two halves of the hysteresis loop from their 
average value at a given applied field. In view of the smallness of the area of the anti-ferromagnetic hysteresis loop, 
the magnified loops are better suited for comparing the fit between simulation and the theory to be presented in the 
following sections. 

IV. BOUNDED DISTRIBUTION OF QUENCHED FIELD WITH A < |J| 

If A < I J| and the applied field ha is increased adiabatically from —oo to cx), the spins with both neighbors down 
flip up flrst (ramp-I). Next are those with one neighbor up and one down (ramp-II). The last category of spins to flip 
up are those with both neighbors up (ramp-III). The ramps are separated by plateau-I and plateau-II. The width of 
the plateaus along the applied field axis decreases with increasing A and goes to zero as A — )• | J|. Magnetization on 
ramp-I was determined in reference [isj by exploiting a similarity between this problem and the problem of random 
sequential adsorption(RSA) (39j . The rate equations of the RSA problem were used to determine m(ha) on ramp-I, 
but they could not determine m{ha) on ramp-II and ramp-III. A different approach was introduced in [l9l. [20| which 
determined magnetization on all three ramps for A < |J|. We recall this approach briefly because it serves as the 
starting point for analyzing magnetization curves for A > \ J\ as well as for a Gaussian distribution of random-fields. 
The following subsections contain the main results obtained in references [19, 20] with some reworking of notation 
and formalism. 



A. Ramp-I 

The analytical results for the magnetization on ramp-I are conveniently expressed in terms of three quantities po(^a), 
Pi{ha), and ^2(^0)- These are the probabilities that a spin which has quenched field hi and which has respectively 
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zero, one, or two nearest neighbors up can flip up at applied field ha- Thus, 



Pn{ha) = / <l>{h,)dh, {n = 0,1,2) (7) 

J-2(l-n)|,/|-/i„ 

As the applied field increases adiabatically from ha = — oo to ha = oo, spins with both neighbors down begin to 
flip up at ha = — 2| J| — A and continue to flip up till ha = J| + A at which stage ramp-I is completed and there 
are no down spins whose neighbors are also down. The fraction of up spins on ramp-I at an arbitrary applied field ha 
is given by. 



P/(;i,) = -[l-e-2^»("")], (8) 

The central object in the calculation of P^(ha) is the probability (per site) of finding a pair of adjacent down spins 
on ramp-I at applied field ha- We denote this object by the symbol P^^{ha)- It is calculated as follows. Imagine 
coloring all sites with hi + 2\J\ + ha > black, and all sites with hi + 2\J\ + ha < white. Consider two adjacent 
down spins A and B shown in Figure (7). The sites A and B can be both white, both black, or mixed. Given that 
A is down, it is clear that the state of B can only be influenced by the evolution of the system to the right of B. 
Similarly, given that B is down, the state of A can only be influenced by the evolution of the system to the left of A. 
Wc shall refer to this as the principle of conditional independence It requires 



Pl^iha)^PiAi\Bi)PiBi\Ai) (9) 

where P{A I \B J.) is the probability that spin at site A is down given that spin at B is down, and P{B | \A l) is 
the probability that B is down given that A is down. We take up the calculation of P{B I \A 1). If B is a white site, 
P{B l\A l) = 1 because white sites have not been relaxed from their initial state. If B is a black site and the site to 
the right of B is a white site then P{B I \A l) = 0. In general P{B l \A |) depends on the length of the string of 
black sites to the right of B. Suppose B is a black site, and there arc {n — 1) additional black sites to the right of B. 
In this case, the probability that B is down satisfies the following recursion relation. 



p^^_p-2^(l__)p«-l (10) 

The rationale for the above recursion relation is as follows. Let the black site farthest from B on the right be labeled 
as the n-th site. Any of the n sites could flip flrst. The probability that the n-th site flips first is therefore equal to i. 
If this happens, {n-l)-th site is prevented from flipping up on ramp-I. The probability that B is down is now reduced 
to the probability that the end point of a chain of {n-2) black sites is down i.e. Pg~^ - This accounts for the first term 
in equation (10). The probability that n-th site is not the first site to fiip up is equal to (1 — ^). Given this situation, 
the probability that B is down is equal to the probability that the end of a string of (n-1) black sites is down. This 
accounts for the second term in equation (10). We can rewrite the recursion relation (10) as 



(p^-pr') = -^[^r'-^r'] (11) 

It has the solution, 



^-^ ml 

m=0 

Summing over various possible values of n with appropriate weight, we get 



P{Bi\Ai) = J2J2 



-Po(l-Po) 



71—0 m— 

oo 



E^(i-^o)EK-E 



i-poY 



(13) 
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In the above array of equations, po stands for po(ha)- Thus, 

Let be the probabiUty per site of finding a down spin and P^^ the probabihty per site of finding a down spin 
which is fohowed by an up spin. Clearly, 

Pi = Pu + Pit = 1 - (15) 

Keeping in mind that on ramp-I an up spin must be preceded (as well as followed) by a down spin, we get P^^ ~ P^ . 
Thus, 

-ft = 1 - Pu ^ Pt 

or, 

The magnetization on ramp-I is given by 

m^{h) = 2P/(/i) - 1 = -e-^PoC^.) (17) 

The exact value of the magnetization on plateau-I is equal to — ^ which is approximately equal to -.135. 

B. Plateau-I 

Plateau-I contains down spins in singlets and doublets punctuated by up spins. Each spin in a doublet has one 
neighbor up and one down. Therefore the net field on it is simply the sum of the random field hi on its site and 
the applied field ha- It turns up when hi + ha > 0. The random field lies in the range —A < hi < +A. Therefore 
an applied field smaller than —A is sufficiently negative to pin down all doublets. This accounts for the range 
(— 2| J| + A) < ha < —A where the magnetization shows a plateau. In each doublet, the spin with the larger quenched 
field hi flips up on ramp-II when hi + ha > 0. The spin with the smaller quenched field then becomes a singlet 
which does not flip up before ramp-III. Thus, in order to find the form of ramp-II, we need to find the a posteriori 
distribution of quenched random fields on the doublets. 

Consider a doublet on plateau-I as shown in Figure (8). The doublet sites are denoted as 1 and 2, and the quenched 
random fields on these sites are hi and h2. The a posteriori probability distributions of hi and /12 will be identical by 
symmetry. These distributions (j){hi) and (j}{h2) are determined by the relaxation process on ramp-I. If the doublet 
survives up to plateau-I it must exist on ramp-I. Consider ramp-I at an applied field ha ■ Given that site-2 is down at 
this point, the probability that site-1 is also down is equal to e~P°'-'*^. Site-1 may be down because (i) it is a white 
site i.e. /ii -|- 2| J| -I- /la < and therefore it could not turn up even if both its nearest neighbors were down, or (ii) 
it is a black site but blocked from turning up by its neighbor that has turned up before it. The probability that it 
is a white site is equal to 1 — Po(^a)- Therefore the probability that it is a black site but down at ha is equal to 
f,-po{h^) _ {l-pQ{ha)]. This means, 

Proh{li\2i-hi + 2\J\ + ha>Q)^ / 4>{hi)dhi^ e-f"(''»)-{l-po(/ia)} • (18) 

J-2\J\-ha '- 

4>{hi) is obtained by taking the derivative of the above expression. We get, 

4>{hi)dhi = [1 - e-P«(-''i-2|J|)]^(^^)^;j^ (19) 

Assuming hi > h2, site-1 would turn up on ramp-II when hi + ha ~ 0. The density of sites on ramp-II at this value 
of the applied field is given by. 
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4>{^ha)dha = [1 - e-^'°(''"-2|^l)]0(-/la)d/ia = [1 - e-P'^^-^^i~ha)dha (20) 

We now address an issue which is crucial for determining ramp-II correctly. This concerns two adjacent doublets as 
shown in Figure (9). What is the probability per site of observing this object on ramp-I? A doublet has an important 
property. It separates the lattice into two parts (one on each side of the doublet) which have evolved uninfluenced by 
each other. Thus, we can separate Figure (9) into three parts as enclosed in the dashed boxes. Evolution inside each 
box has remained shielded from the outside. The evolution in the middle box requires that site 3 flips up before site 
2 or site 4. The probability for this event is equal to i. Given that site 2 is down, the probability that site 1 is down 
is equal to e~'P"^'^'^K Similarly, given that site 4 is down the probability that site 5 is down is equal to e"^"''''"). Thus 
the probability per site of observing two adjacent doublets on ramp-I at an applied field ha is equal to ie~^^°'-''°'. 
Note that it is quite different from the square of the probability of finding a single doublet! Let hi, h2, /13, /i4 and /15 
denote the quenched fields at sites 1, 2, 3, 4, and 5 respectively. We are interested in the case /12 > hi, and /14 > /15, 
and ask what are a posteriori distribution of fields /ii, /12, • • • , on ramp-II. Let the probability that hi{i = I, . . . ,5) 
lies in the range [—ha — dha, —ha] on ramp-II be denoted by pi{—ha)dha. The distributions of hi and are given by, 

Pli-ha)^ ^ ^ H~ha) (21) 

Pd~ha)^ ^ ^ 0(-M (22) 

Equation (21) is obtained as follows. If sites 1 and 2 are down on plateau-I. they must have been down all along 
ramp-I. Given that site 2 is down on ramp-I, the probability that site-1 is down is equal to e~P°^^'^\ This accounts 
for the denominator. The numerator gives the a posteriori probability distribution of the quenched field at site 1 on 
the lines of the preceding discussion with reference to figure (8). Equation (22) is written similarly. Note that the 
distributions pi{hi) and p^ihi) are each normalized to unity. Next we turn to the distributions of /i2, /is, and /14. 
As stated previously, /12 > hi and h^ > so we are looking at the case where ft.2 and h^ flip up on ramp-II. We 
may assume without loss of generality that > ft,2- Let flips up on ramp-II at applied field ha i.e. + ha — 0. 
The probability for this is equal to c/)(— /ia)pi(/ia){l ~ Pi{ha)}', the three multiplicative factors giving respectively the 
probability that /i4 — —ha, site 3 is up at ha, and site 2 is down at ha- Thus, 

Pi{-ha) = 6pi{ha){l - Pl{ha)}cl){-ha) (23) 

The factor of 6 on the rhs arises as follows. At ha on ramp-II site 2 will be down and sites 3 and 4 will be up at 
ha given that at an earlier field ha — 2\J\ on ramp-I site 3 had flipped up before 2 and 4 with probability i. The site 
flipping up on ramp-II may be to the left of the central site or to its right. This gives an additional factor of 2. 

The distribution of ft,2 is obtained similarly. We find, 

P2{-ha) = 3pl{ha)(b{-ha) (24) 

C. Ramp-II 

Ramp-II is determined by the combination of two terms. The dominant term is the increase in magnetization due 
to the decrease in the number of doublets. When a doublet disappears, it adds an extra up spin in the system which 
increases the magnetization. Occasionally, a disappearing doublet creates a string of three up spins. A triplet of up 
spins is unstable on ramp-II if A < | J| and therefore the central spin of the triplet flips down as soon as the triplet 
is created. This decreases the magnetization. In the following, we calculate the above two terms separately. Refer to 
figure (8) for calculating the first term. Let us assume hi < /i2- The a posteriori distribution of fields hi and /i2 in 
figure (8) are the same as pi{hi) in figure (9). The probability that the doublet disappears at h2 + ha — on ramp-II 
is given by 

/oo /-ho 
Pi{h2)dh2 / " Pi{hi)dhi (25) 
-ha J —OQ 
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The factor e"^^''^''"-' is the probabihty per site of finding a doublet at ha before any of them have been relaxed. 
The factor 2 takes care of the fact that either hi or /i2 may be the larger field although the expression is written on 
the assumption that site 2 flips up first. When a doublet disappears, a pair of adjacent up spins is created. This is 
the reason for the choice of the subscript on . The superscript indicates that the probability refers to ramp-II. We 
obtain, 



h2 — oo 



h2— — ha 



(26) 



For A < |J|, po{ha) = 1 if pi(/ia) > 0. Therefore, 



pii - ± 



(27) 



We now calculate the fraction of (unstable) up triplets on ramp-II. Refer to figure (9) with the assumption that 
/i2 < hi. An up triplet forms when h2 + ha = 0. The cumulative fraction of up triplets is given by, 



1 



P//^ = -e-2Po(/>») / p^{h2)dh2 pi{hi)dhi pi{hi)dhi P5{h5)dh5 

"J J -h„ J -oo Jh2 J-oo 



(28) 



The two factors before the integrals give the probability per site of finding the object shown in Figure (9). These 
take into account the condition that /13 > /14. The next two integrals give the probability that h2 + ha = and 
hi < h2. When ft,2 is in the range —ha and —ha — 5ha, /14 can be anywhere in the range /12 to 00. Let p4(hi) be the 
density of /14 in this range. Clearly, 



P2ih2 



Pi{hi)dhi, or p4,{hi) 



dp2{hi) 
dhi 



(29) 



The integrals in equation (28) can be evaluated exactly for a uniform distribution. We get a non-zero contribution 
only \i pi{ha) > 0. If A < | J| and pi(ha) > then we must necessarily have po{ha) = 1. Thus we get, 



^ + ^ - 6 |l + l\pi{ha) + Spjiha) + PU^a) " + l} + 



6 (^1 + - 6pi{ha) - 3 (^1 + pjiha) + 2p?(/ia)| e-P'^'''^ 1^ + 3pi(/ia)| e-^PiC'^) 



Putting the various terms together, the probability that a randomly chosen spin on the lattice is up on ramp-II is 
j'iven by 



(ha) = ^[l - e-'] + Piiiha) - P^-^{ha) 



(30) 



The magnetization on ramp-II is given by 



m'\ha)^2PJ^'{ha)-l 



(31) 



The magnetization on platcau-II is equal to. 



m 



II 



27 7 _i 8 _2 

e ^ e 

30 6 3 



.109 (approximately) 



(32) 
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D. Plateau-II 



Each down spin on platcau-II is a singlet. However, there are three different classes of singlets: the singlets formed 
on ramp-I; singlets formed on ramp-II by a vanishing doublet; and finally the singlets formed on ramp-II by the 
unstable central spin of an up triplet flipping down. The a posteriori distribution of random field is different for each 
class. Let p^J , pl^ , and p^^ denote the density of random fields in the three cases and P^^ , P^^ , and P/^ be the 
probability of finding the corresponding singlet at applied field haOn plateau-II. 



p'J{hi)dh.i 



b, c 



(33) 



It is useful to think of singlets in each class as being black or white on the ramp on which they are created. Suppose 
a singlet is created on ramp-I at applied field h'^. It is created black ii hi + 2\J\ + h[^ > 0, and white ii hi + 2\J\ + h'^ < 
where hi is the quenched random field at the singlet site. If a singlet is black at applied field h'^ then it is also black 
at fields greater than h'^. A black singlet on ramp-I at /i^ will turn up on ramp-III at ha = h'^ + 4| J|. Therefore if we 
know the fraction of black singlets on ramp-I we can calculate how they are destroyed on ramp-III. 

The fraction of black singlets created on ramp-I is given by, 



Pl\K) = p,{h',) - - [l - e-'P^^'^'^A - 2e-f°(''^) [e-^«(''^) - {1 - p,{K)} 



(34) 



The explanation of the above equation is as follows. Imagine ordering the sites of the lattice in order of decreasing 
quenched field on the site. When all sites with hi + 2\J\ + h'^ > have been relaxed, the fraction of the relaxed sites 
is equal to po{h'^) (the black sites). This fraction is made of the up sites (the second term on the right), black doublet 
sites (the last term), and the black singlets. Hence the equation for P(^{h'^). The last term is written as follows. In 
each doublet, there are two sites from which we can choose one. This accounts for the factor 2. The quantity in the 
square bracket gives the probability that the chosen site is black, and e^^"'-''"-' is the probability that the other site 
can have any allowed value of the quenched field. 

It is instructive to derive equation (34) by an alternate and more direct method as well. We note that a pair 
of down spins on the chain has to be followed by a down spin or an up spin, i.e. P^ = Pm + Pii-\-- Similarly, 



P^ 



p 



Pu- Thus, 



PniiK) - PmiK) = PuiK) PmiK) = e-'^"^"-^ - {l - Poih',)}e-'P"^'^'^^ = poih',)e-^PoiK) 



(35) 



We can obtain the fraction of singlets P^i^ from P^n and Pn^- by calculating the probability that the down spin 
at either end flips up under a po process. 



PntiK) - 2 / p^uiK)<l^i-m - K)dh: = 2 



-2po 



-2po('i'o: 



Po(/i;)e-^P«(''") (36) 



The above expression gives the total fraction of singlets on ramp-I which include black {hi -f 2| J| + h'^>Q) as well 
as white [hi + 2\J\ + h'^ < 0) singlets. The fraction of white singlets is given by. 



^m*^(/*l) = {i-po(/il)}[i-e-^'«(''»)]2 



(37) 



The first factor on the rhs gives the probability that the central site is white. Given that the central site is white, 
the probability that its neighbor is up is equal to 1 — e"^"^''"-'. The probability that both neighbors are up is equal 
to [1 - e-P«(''^)]2. Thus, 



P^'C^HK) = PmiK) - Pm'%K) = PoiK) - 2 [l - e-^^^e^")] - 2e-P<'^<^ [e-P"^'^'^^ - {1 - poih'J} 
The fraction of black singlets on ramp-II {hi + /i" > 0) generated by vanishing doublets is given by. 



(38) 



Pi'iK) 



-piiK) 



-{1-pM} 



(39) 
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The above equation is easily understood. It is the probabiUty that both sites of the doublet are black. If both 
sites of the doublet are black, the one with higher random field must flip up on ramp-II leaving us with a singlet on 
plateau-II which is black. 

The fraction of black singlets created by unstable triplets requires the calculation of triplets. We have calculated 
the fraction of triplets as they are formed on ramp-II. What wc need now is a similar but different calculation. The 
point can be understood with a reference to Figure (9). Recall that /13 > /14 > /i2- On ramp-II, we needed the fraction 
of triplets with /12 + /i" > 0, because the formation of triplets is controlled by this threshold. The restoration of the 
triplets on ramp-III is controlled by the condition — 2\J\ + ha > 0. Keeping in mind that we want hi < /12, and 

< h4, the probability that /13 — 2| J| -I- /la > is given by, 

/oo ph^i fhi phi 

(j){hs)dh3 / p4{h4)dh4 / P5{h5)dhz / (02(^2)^^2 / pi{hi)dhi 
-/la+2|J| J — 00 J —00 J —00 J — 00 

Using po{ha) = 1 and pi{ha) = 1 on platcau-II we get, 



i [1 + I] Ptiha) + ^pUha) +[! + §] e-P-^'^'^ - |p2(/^a)e-^'^("") +p2(/i,)e-P^(''") 



= -2p2('>a)l 



(40) 



E. Ramp-III 

The rise of magnetization on ramp-III is due to singlet sites turning up in increasing field. At the start of ramp-III 
there are three categories of singlets present on plateau-II, and we have classified each of them conveniently into black 
and white singlets. The fraction of singlets on plateau-II that turn up at ha on ramp-III is given by the fraction of 
black singlets in each of the three categories: 



P''\ha) = Pi" (ha) + Pi" (ha) + P^" (ha) (41) 

However the calculation of magnetization on ramp-III turns out to be a bit more complicated. There is a new twist. 
Frequently when an original singlet site on plateau-II turns up on ramp-III its nearest neighbor turns down. We call 
this the creation of a new singlet on ramp-III. The newly created singlet site would turn up at a larger applied field on 
ramp-III. We shall call this event the destruction of the newly created singlet. We have to calculate the newly created 
singlets and their destruction before the magnetization on ramp-III may be obtained. It should be noted that this 
means that some sites flip three times in the course of a monotonic increase of applied field from — oo to oo. However, 
no site flips more than three times. 

When does a vanishing singlet on ramp-III create a new singlet on an adjacent site? Consider the singlet at site 3 
in flgure 10. Suppose site 3 flips up at ha, i.e. /13 — 2\J\ + ha = 0. Now the net field on site 2 is equal to /12 + ha- 
This is necessarily positive because /12 + 2| J| — /13 > if A < \ J\. Thus site 2 would stay up after site 3 flips up. 
Consider site 4. After site 3 has turned up the net field at site 4 is equal to /i4 — 2| J| -I- /iq. Thus site 4 will turn 
down if /i4 < hj. Site 3 will stay up even if site 4 turns down because h^ + ha > 0. These considerations can be 
put in the form of two guiding rules. When a singlet turns up on ramp-III, (i) its nearest neighbor stays up if the 
next nearest neighbor is down, and A < | J|, (ii) its nearest neighbor turns down if it has less quenched field than the 
singlet and the next nearest neighbor is up. Detailed considerations show that only the singlets created on ramp-I 
fall under the purview of these rules. Therefore we focus on the singlets present on plateau-I. Specifically we focus on 
the configurations shown in figure (11) and figure (12). In each of these figures a new singlet is created on site 3 when 
the singlet on site 2 is destroyed. The two figures make different contributions because of the role played by the next 
nearest neighbor of the singlet on the side of the newly created singlet. 

The a posteriori distribution of the quenched field at site 2 is given by 

kh2) - [1 - e-P^(-''^+2|''l)] 0(/i2) (42) 

The contributions of figure (11) to the fraction of newly created singlets when all sites with hi ~ 2\J\ + ha > have 
been relaxed on ramp-III is given by. 
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4>{h2)dh2 / 0(/l3)d/l3 / 0(/l4)d/i4 

-ha+'2\J\ J -oo Jh2 

^ - 2p2{l-p2) - \pI - 2(1 -P2 +pl)e-'''' + - 2P2)e''^'^ [p2 = P2iha)]. 



(43) 



Similarly the contribution of figure (12) is, 



/oo _ i<h2 rhs i-hi _ 

4>{h2)dh2 / 0(/l3)d/j3 / 0(/i4)dM / ~^{h5)dh5 
-/la+2|J| J -OO J -OO J -OO 

= - + 3e-i) + Q + 5e-i) P2 - + 2e-i) + ^ (l + e-^) 
-^P^ + { + 3e-i) - (1 + 2e-i) + - Ip^} e"''^ - e'^^^ 

(44) 

In calculating the total fraction of newly created singlets, we have to multiply the contribution of figure (12) by 2 
because an equal contribution is made by a configuration in which the doublet is to the left of the vanishing singlet. 

The destruction of newly created singlets on ramp-III can be analyzed in a similar manner as their creation. The 
probability that site 3 flips for the third time in figure (11) is given by. 



Py'{ha) = 2 Cf>{h3)dh3 Hh2)dh2 ^{h4)dhi 

-/l„+2|J| J/l3 Jh3 

^[l^e-'P^]~2p2e-P^+P2{l-p2) + lpl 



(45) 



Similarly, the contribution to destruction of newly created singlet in figure (12) and another similar figure in which 
the doublet is to the left of the singlet is given by. 
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P^'\ha) = 2 Hh3)dh3 (j){h2)dh2 ^ihi)dh4 (t>ih5)dh5 

J-ha+2\,J\ Jh3 J-oo J-oo 

= 2e-i - (1 + 4e-i) P2 + Q + Se-^^ pI-[i + ^e-^^ pi + \pi 
- { (1 + 2e-i) - 2 (1 + e-i) P2+pl) e'^'^ + e^^^' 

(46) 

Putting all the terms together, the probability that a randomly chosen site on ramp-III is up is given by, 

P/"(/la) - P^\ha) + Pi" (ha) + Pl'\ha) + P'J' {ha) - Pj' ' (ha) - P^'^ha) + Pj" [ha) + P/"(/ia) (47) 

The magnetization in increasing field ha is given by m(ha) = 2P^^^ [ha) — 1. The magnetization m/j on the 
return trajectory in decreasing field may be obtained by symmetry mniha) = —m{—ha)- These results has been 
superimposed on the corresponding simulation data in figure (1) and figure (2). The simulation results for the return 
hysteresis loop were obtained independently without using the symmetry. 
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V. UNBOUNDED DISTRIBUTION OF QUENCHED FIELD 



It is convenient to introduce the following nomenclature. We say that a site flips up under a po-process if none of 
its two nearest neighbors are up when it flips up. It is said to flip up under a pi-process if it flips up when one of its 
neighbors is up and the other is down. Similarly, a site flipping up under a p2-process has both its neighbors up at 
the time it flips up. The simplifying feature of the analysis for A < |J| is that a p„-process (n = 1,2) can not take 
place anywhere on the chain unless all Pn-i processes have been exhausted. In other words a new ramp can not begin 
before the previous ramp is completed. This feature is lost if A > | J| (or if the distribution is Gaussian) because a 
Pi or a p2 process can occur on the chain even if po-processes have not been exhausted (i.e. there remain strings of 
three down spins on the chain). The fact that po, pi, and p2 processes can run concurrently makes the calculation of 
a posteriori distribution of quenched fields at down sites more complicated than encountered in the preceding section. 

Our first task is to calculate the probability of occurrence of a doublet in the chain. A doublet is a pair of adjacent 
down spins that have remained down in a monotonically increasing field from — oo to ha- The significance of a doublet 
lies in its screening property. It separates the chain into two parts that have evolved independently of each other. The 
doublet also provides a natural length in the analysis of the chain. This is not a fixed length but rather a variable 
length of a segment of chain that is free of doublets and lies between two doublets at an applied field ha- The history 
of evolution of this segment may be analyzed independently of the rest of the chain. Let us focus on one end of 
such a segment. We label the sites 0, 1, 2, 3, . . . , n starting from the left end of the segment. Our immediate object 
is to calculate the probability (per site of the chain) that site-0 and site-1 form a doublet at ha- This is equal to 
P^j^(l i |0 i;ha) where Pii{l i |0 i;ha) is the conditional probability that site-1 is down given that site-0 is down. 

We proved in section IV that ^4,4.(1 i |0 I; ha) ~ e"^"'-''''' if only pg-processes are allowed. This is an exact result on 
ramp-I for A < | J| and we take it as the leading term of the exact result for A > | J| or for an unbounded distribution 
such as the Gaussian distribution. The main effect of A < \J\ is that if one of two adjacent sites say site-1 and site-2 
{hi + 2\J\ + ha > and /12 + 2\J\ + ha > 0) flips up under a po-process then it prevents the other from doing the 
same. As long as \hi — ft.2| < 2\J\ the same effect will be obtained at site-1 and site-2 for other distributions of the 
quenched field. In any particular realization of the quenched field distribution, the occurrence of a large connected 
cluster of sites with hi+i — hi> 2\J\ is rare. Therefore the case A < | J| serves as a good starting point for the exact 
result. Thus, 



Pu(l i |0 ha) = e-P^C*") { leading term } (48) 

Our approach is to add correction terms to the leading term to make it an exact result. The correction terms are 
functions of pi{ha)- No correction to the leading term is required ii pi{ha) = 0. For pi{ha) > 0, the corrections may 
be divided into two categories. The first category is the one in which a doublet created by a po-process is subsequently 
destroyed by a pi-process. This is similar to the fate of doublets on ramp-II for A < | J| and we get, 



P;4(l ; |0 i; ha) e-P«(''") - L-Pi^'*") - {1 - PiiK)}] { first order correction } (49) 

The second category of correction involves events where a pi-process pre-empts a po-process. These events are 
a signature of A > |J| and do not exist if A < |J|. Such events are always possible if the quenched field has a 
Gaussian distribution. We illustrate this by a simple example. Suppose for a particular realization of the distribution 
of quenched fields ft-i, /i2, h^^ /14 . . ., site-1, site-2, site-3 are black, site-4 is white, and h^ > h2 > hi- If A < | J|, site-3 
will flip up first, block site-2 from flipping next, and therefore site-1 will flip up last. This will result in the absence of 
doublet at site-0 and site-1. Now consider A > \ J\ so that it is possible to have /i2 > hi + 2\J\ but hi + ha < 0. In this 
scenario site-3 will flip up first, then sitc-2, and site-1 will remain down giving us a doublet at sitc-0 and sitc-1 that 
was not allowed if A < | J|. The probability of this new doublet created by a pi-process pre-cmp tying a po-process is 
given by, 



/—ha pOO 
(j){hi)dhi / ^{h2)dh2 (50) 

-2\J\~ha Jhi+2\J\ 

Here 0(/i2) is given by equation (19) of section IV. This is because equation (19) of section IV gives the a posteriori 
distribution of the quenched field on a black site next to a site that has fiipped up by a po-process. In our example 
site-2 is a black site next to site-3 that has flipped up by a po-process. T^{ha) term is the leading term in the category 
of correction terms that arise because a pi-process has pre-empted a po-process. 



14 



Pu{li\Oi;ha) 



e-P'^'^'^ ^ {1 - pi{ha)}] +niha) 



{ second order correction } 



(51) 



The next correction comes from a cluster of adjacent spins that flip up as follows. Site-4 flips up under a po- 
process (ft.4 + 2\J\ + ha > 0), site-3 flips next under a pi-process {h^ > /12 + 2| J|), site-2 also flips under a pi-process 
(/i2 > /ii + 2| J|) after site-3 has flipped up, and site-1 remains down {hi + ha < 0, hi+2\J\+ha > 0). The contribution 
from this event is, 



-ha roc pQO 

Ti{ha) = - / (p{hi)dhi / ^{h2)dh2 / ^(/l3)d/l3 

-2\J\-ha Jhi+2\.J\ Jh2+2\J\ 



(52) 



Notice that T^{ha) is positive and Ti{ha) is negative. T^{ha) is positive because it produces a doublet (at site-0 
and site-1) where it did not exist under the po-process alone. Ti{ha) is negative for the following reason. r4-process 
and po-process is mutually exclusive and since both contribute to a doublet then they must be added separately 
avoiding double counting. Thus the leading term in equation (48) under po-process alone is an overestimate and has 
to be reduced by an amount T/^{ha)- However, r4-process also gives rise to a doublet. Hence an amount equal to 
the one subtracted from the leading term has to be added to it resulting in zero correction to the leading term. Now 
consider the first correction shown in equation (49). This is an underestimate because it comes from the destruction 
of overestimated doublets in the leading term. It can be corrected by adding T^lha). The correct result at the present 
level of accuracy is. 



Pu(l ; |0 i; ha) = e-P"^''-^ - [e^PiC'") - {1 - pi{ha)}\ + niha) - Ti{ha) { third order correction } (53) 

Continuing in this vein we get the following exact result, 



-Pliha) 



{l~Pi{ha)}] - ^(-l)"r„(/i,) { exact result } 



n=3 



where 



Tniha) = / (j){hi)dhi 

-2\J\-ha 



'n-1 

n 

.m— 3 



4){hrn-l)dh 



m — 1 



h^^2+2\,J\ 



(f>{hn-l)dhn- 



(54) 



(55) 



/l„_2+2|J 



The probability per site of finding a doublet on the chain is equal to, 



Puiha) = I |0 ha)] 



(56) 



Note that although the (unconditional) probability of a doublet is the product of two mutually conditional prob- 
abilities, our notation for this is Pii{ha) (without the square sign). In general it is not possible to do the integrals 
exactly to get an analytic expressions for T„(ft,a) in a closed form. These have to be evaluated numerically. However 
Tn decreases exponentially with increasing n. In our numerical work we included terms up to n < 4 in the calculation 
of the conditional probability Pii{l i |0 j,: ha), and used the square of this quantity to calculate the probability per 
site of doublets. This gives an excellent fit with the simulation data as shown in section VI. 

Although a doublet at an applied field ha is the most basic object in our calculations but we need to calculate 
the probability of several other objects before we can calculate the magnetization curve. In order to calculate the 
magnetization curve we need to know the probability P-^{ha) that a randomly chosen site is up at ha- On the lower 
half of the hysteresis loop, it is more convenient to focus on the complementary probability [l-P-|-(/ia)] that a randomly 
chosen site is down. The randomly chosen site can have both its neighbors down, or both of them up, or one up and 
one down. If both neighbors of a down site are down, the site in question is necessarily a white site. The probability 
(per site) that a site is white is equal to 1 —poiha)- Given a white site, the probability that both its neighbors are 
down is given by. 



Pmiha) = [1 - Po{ha)]Pii{ha) 



(57) 
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Next, consider a down spin with one neighbor up and one down. In both cases there is a doublet, and the up spin 
is either to the left of the doublet or to its right. We have calculated the probability of a doublet as well as a triplet. 
From these we can obtain the probability of a doublet followed by an up spin using the equation 



Pmiha) ^ Puiha) - Pmiha) (58) 

Naturally Pufiha) = Puii^a) by symmetry. 
Substituting from equations (57), we get 

Pmif^a) = Pniiha) = P0iha)Puiha) (59) 

Note that the up spin in the above object could have flipped up under a pa or a pi process. We have checked the 
above equations against numerical simulations shown in section VIII. 



VI. SINGLETS 



Now we calculate the probability P^i-f{ha) that a randomly chosen site is down and both its neighbors are up. Let 
the three consecutive sites that form this singlet be labeled 1, 2, 3 and hi, /12, and /13 denote the respective quenched 
fields. Without loss of generality we can assume that just before this singlet is created site-1 was up and sites 2 and 
3 were down as shown in figure (14), i.e. site-1 flips up before site-3. We have to know if site-1 has flipped up under 
a po-process or a pi-process. We consider both these possibilities for site-1 as well as site-3. The case when site-2 has 
never flipped in the course of applied field changing from —00 to ha is somewhat simpler to analyze. Even in this case 
the probability P-^i-f(ha) depends on whether (i) both neighbors flipped up under a pg-process, or (ii)both neighbors 
flipped up under a f>i-process, or (iii) one neighbor flipped up under po-process and the other under a pi-process. In 
the case (iii) it is also important whether the neighbor that flipped up first flipped under a po-process or a pi-process. 
Indeed we require the following objects before we can calculate P-fi^{ha) in the simpler case mentioned above. 

• P-ti^A The fraction of singlets at applied field ha when site-1 is up by a po or a pi-process and site-3 flips up by 
a po-pi'ocess. 

• -PfJ-tB Fraction of P-\-i-\-a that are black at creation. 

• ^tJ-tc Fraction of P-\_i^a that are white at creation. 

• PtitD The fraction of singlets at applied field ha when site-1 is up by a po or a pi-process and site-3 flips up 
next by a pi-process. 

• ^^tJ-t^; Fraction of P^i^d that are black at creation. 

• P-\i-tF Fraction of P^i^d that are white at creation {P^i^d — -Ptit^^)- 

• P-fltG The fraction of singlets at applied field ha when site-1 and site-3 flip up by a po-process and site-2 is 
white. 

• P'ti'tH The fraction of singlets at applied field ha when site-1 flips up by a pi-process, site-3 by a pg-process, 
and sitc-2 is white. Site-1 flips up before site-3. 

• P-tl"ti The fraction of singlets P-fi^c that are white at their creation but are black at ha- 

• ^tJ-t-' The fraction of singlets at applied field ha when site-1 and site-3 flip up by a pi-proccss and site-2 is 
white. 

• P-\i^K The fraction of singlets at applied field ha when site-1 flips up by a pi-process, site-3 by a pg-process, 
and site-2 is white. Site-3 flips up before site-1. 

• P-\l^L The fraction of singlets P-fi^fp that are white at their creation but are black at ha- 

• ^tttM Destruction of singlets associated with P-\i-[b- By destruction we mean the disappearance of a singlet 
due to the down spin flipping up under a p2-process. 

• PtttN Destruction of singlets associated with P^i^e- 
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• -Fttto Destruction of singlets associated with P-|-4.-|-c- 

• ^ttt^^ Destruction of singlets associated with Pfifp- 

Up to this point we have focused on singlet sites that have never flipped starting from the saturated state at 
ha = — oo. In other words we have considered configurations comprising site-1 up, site-2 down, site-3 up where site-2 
has never flipped up in increasing field from — oo to ha- However, the anti-ferromagnetic dynamics also allows a singlet 
with sitc-2 having flipped twice i.e. site-2 can flip up and flip down again in the course of monotonically increasing 
applied field. In order to take into account the analysis of this second category of singlets, we need to calculate the 
following objects that are defined with respect to figures (15), (16), and (17). 

• -^tJ-tQ This object refers to figure (15) with the proviso that /12 > hi, /14 > /15, /14 > /12 (/14 < /12 will give an 
equal contribution). Site-3 flips first, sitc-4 fiips next, site-2 flips after site-4 causing site-3 to flip down. P^i^q 
refers to the fraction of singlets created in this way. 

• ^tJ-tfl This object refers to figure (16) with the proviso that /12 > h^, /14 > /13, /14 > /12 (hi < /i2 will give an 
equal contribution). Site-3 fiips first, site-4 fiips next, site-2 flips after site-4 causing site-3 to flip down. P\i-[r 
refers to the fraction of singlets created in this way. 

• ^t-l-t'S This object refers to figure (17) with the proviso that /12 < /13, /14 > /13. Site-3 fiips first, site-2 fiips next, 
site-4 flips after sitc-2 causing sitc-3 to flip down. Pfifs refers to the fraction of singlets created in this way. 

• Ptl1-T This object also refers to figure (17) with the proviso that /12 < /13, /14 > /13. Sitc-3 flips flrst, site-4 flips 
next, sitc-2 flips after site-4 causing site-3 to flip down. P^i-\-t refers to the fraction of singlets created in this 
way. 

• ^tttc^ This object refers to the destruction of singlets associated with P^i^q- 

• ^ttt^^ This object refers to the destruction of singlets associated with P^i^b- 

• Ptttw This object refers to the destruction of singlets associated with Pfifs- 

• ^tttA' This object refers to the destruction of singlets associated with P-|-4,-|-t- 

A. P-f^-fA 

In an increasing applied field the objects associated with Pfi^A are created from objects associated with P^f and 
P-tii when the down site at one end of these objects flips up under a po-process. Suppose site-3 flips up under a 
Po-process at h' (—00 <h'< ha)- Then we have /13 -|- 2\J\ + h' ~ Q with probability 4>{—2\J\ — h'). Thus, 

PmA{ha) - / ° [Pm{h') + Pni{h'm-h' - 2\J\)dh' - / " 2po{h')P^^{h')(^{-h' - 2\J\)dh' (60) 

J —oo J —oo 

where we have used equation (59). As a check we note that equation (59) is recovered by differentiating equation 
(60) with respect to ha- 



B. P-Xl^B 

This quantity is given by the equation. 



PmB{ha)=2 j \p^^{2i\ii:h')-{l-po{h')}]P^^{^i\?,i:h')(^{-h' -2\J\)dh' (61) 

The explanation of the above equation is as follows. The integrand comprises three factors. The last factor is the 
prob that site-3 flips up at h' under a po-process. The other factors take into account that site-2 and site-4 (the right 
neighbor of site-3) arc down at h' and site-2 is black: 4){—h' — 2\J\)dh' is the probability that site-3 flips up at h'\ 
Pll{2 i |3 h') - {l-pa{h')) is the probability that site-2 is down and black (/i2 + 2\J\ + h' > 0); P^^ii i |3 i; h') is 
the probability that site-4 is down. 
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C. Ptitc 

This is simply equal to P^i^a — PtlfB but the notation is useful in the foUowing analysis. 



Pmc{ha) = PlitA{ha) - PntBiK) 



(62) 



D. R 



The singlets in this category are generated by doublets that are bordered by up spins at both ends. Thus site- 
1 is up, site-2 and site-3 are down, and site-4 is up. Site-2 and site-3 are on equal footing and therefore the a 
posteriori distribution of quenched fields 4'{h2) and (j){h^) are identical. We can assume without loss of generality that 
/13 > /i2 and multiply the result by a factor 2. Thus we focus on a singlet created at sitc-2 when sitc-3 flips up at h' 
(-00 < h' < ha). We get, 



Pnwiha) = 2 / dh'^i-h')p„ih')Pui2 ; |3 h') 



(63) 

The explanation of the above equation is as follows: <j){—h') is the probability that site-3 flips up at h' + = 0; 
Po{h')Pn{2 l |3 I; h') is the probability that site-2 is down and site-1 is up just before site-3 flips up. The a posteriori 
distribution (/)(/i') may be obtained by differentiating the following equation. 



dhi4>{h:i)= \e-^-^'^'^ -{l-p^{h')}\-Y,{-^rt^{h' -2\J\) 



71 — 3 



where, 



(64) 



Xn{h') = / cty{h^)dh^ 

J~2\J\-h' 

We may rewrite equation (54) as 



'n-l 

n 

_m— 3 



/i,„_2+2|,7| 



/l„_2 + 2|,/| 



0(/l„_l)d/l„- 



(65) 



Pu(l ; |0 h') = e"P°(''') - [e-fiC'') - {I - p^{h')}\ - ^^(-l)" [f„(/.') - f„(/j' - 2\J\) 



n=3 



(66) 



For spins flipping up under a pi-process, only the second and the last term on the right hand side come into play. 



E. P^ 



Here we want the fraction of singlets in section D that arc created black. If h' is the field at which site-3 flips up 
by a pi-process then the cumulative fraction of singlets which are black at creation is given by 



PmEiK) = 2 / dh'^{-h') [Fu(2 ; |3 h') - {1 - p^{h')}] (67) 

J —OQ 

The first factor in the integrand gives the probability that /13 + h' = 0, and the second factor gives the probability 
that sitc-2 is down and black at h'. 



F. P^lfF 



PfifFi^'^a) = PntD(ha) - PufEiK) 



(68) 
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t-l-tG 



If a site is white with probabihty 1 — po(^a) then it is down with probabihty unity. The probabihty that either of 
its neighbor is down under a po-pi'ocess alone is equal to e~P°^'^''\ Thus the probability that either of its neighbor is 
up under a po-pi'ocess alone is equal to 1 — e^P'>'-^°-\ Therefore, 



Pmciha) - {1-P0{ha)}[l ~ e-P°^'^-^f 



(69) 



H. 



PmHiha) = 2{1 - poiha)} / dh'Pu{3 i |2 h')c^{-2\J\ - h') / cp{^h")dh" 



h' 



(70) 



The above equation is understood as follows. The integral over h" takes care of the site flipping up under a pi-process 
at h" . Therefore the density associated with this integral is the a posteriori distribution 0(— ft,"). Site flipping up 
under the po-process flips up at h' where h' > ft". The density associated with this site is P|,4,(3 J, |2 J,; ft')0(— 2| J| — ft.'). 
The first factor accounts for the fact that the site in question is down given that it is next to a down site (the white 
site with probability 1 — po{ha) and the second factor accounts for the fact that although it is down it is on the verge 
of turning up at ft' under a po-process. Finally the factor 2 takes care of an equivalent configuration in which the 
locations of sites flipping up under a pq and a pi process are interchanged. 

I- P-fiti 

The fraction of singlets P^i^c that are white at their creation but black at ha is given by, 



(71) 



J. Pi 



tit J 



Pm-ri^a) = [1 -po{ha)] 



dh^4>{h'i) 



(72) 



The first factor takes into account that the singlet site is white. The second factor gives the probability that both 
neighbors of the down site flip up under a pi-process. 



^'tit^^ '^^y obtained on similar lines as P^\,-\h- We get, 

PmK{ha)=2{l-po{ha)} r°dft'<^(-ft') r dft"Pu(3;|2;;ft")<^(-/i"-2|J|) 



(73) 



The fraction of singlets P^^^f that are white at their creation but black at ft^ is given by. 



(74) 
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M. 



We now turn to the destruction i.e. the disappearance of a singlet due to the down spin flipping up under a 
P2-process. We start with the destruction of singlets associated with P^i^^b- Taking into account the inequalities 
/i2 — 2| J| < /13 < /i2 that require sitc-2 to be down and a black site when site-3 flips up, we get 

PntMiha) = 2 / dh'4>{-h' + 2\J\) / dh"(l){-h" - 2| J|)Pu(3 i |2 i; h") (75) 

J-OD Jh'-4\.J\ 



N. 



We now turn to the destruction of singlets associated with P-]^i^e- These singlets were created under the condition 
h2 + 2\J\ > hs > h2. They will be destroyed at /12 - 2\J\ + h' = with the probability (j){-h' + 2\J\). Thus the 
cumulative fraction of the destroyed singlets at ha is given by, 

l-ha i-h'-2\J\ 

PntN{ha) = 2 dh'4>{-h' + 2\J\) dh"4>{-h") (76) 

J-oD Jh'-4\J\ 



o. 



ttto 



We now consider the destruction of singlets associated with P-|-4.-|-c- Recall that Pfifc are white at creation. A 
fraction P-xi^i of these become black, say at h' i.e. these can flip up at h' under a po-process if they were to have both 
neighbors down. However they have both neighbors up. Therefore they would flip up at applied flcld ha = h' + 4:\J\. 
Thus we get, 

Pmo^PmiC^a-^lJl) (77) 

F*. .P-]--]-]- p 

P^ttP represents the destruction of singlets associated with P-\i-fF- This is similar to the preceding case because 
P-\l-\-F arc also white at creation. Following similar reasoning as used in the preceding section, 

Pmp ^ PntLiha - i\J\) (78) 



Q. 



t-l-tQ 



Refer to flgure (15) and the deflnition of -Pf^tQ. We have hi < /12, < /14. The screening property of doublets 
ensures that the evolution of sites 2,3, and 4 is uninfluenced by sites 1 and 5. Also, /13 > /14 > /12 for site-3 to have 
flipped up before sites 2 and 4 (/13 > /12 > h^ would make an equal contribution). Now site-2 flips up at ha and site-3 
flips down. Therefore, h2 + ha ~ 0, and /13 — 2\J\ + ha < 0. Consequently /13 — /12 < 2\J\, /14 < /13 < /12 + 2| J|, and 
/i2 < /14 < /i2 + 2| J|. The probability that site-3 would flip down when site-2 flips up is given by 

PmQiha) - 2 dh'4>{~h')p^^{i ; |2 h') f dh"4>{-h")p^^{^ ; |4 h") f dh"'4>{^h"' - 2\j\) (79) 

J-00 Jh'-2\.]\ J-h" 

The limits on the integrals were discussed just before the equation. Note that one can use the quenched field at a 
site as a variable of integration or equivalcntly the applied field at which the site in question fiips up. We have written 
the integrals in terms of the applied fields h'" , h'\ h' at which sites 3,4, and 2 flip up respectively ( h'" < h" < h'). 
The explanation of the integrands is as follows. The integrand in the last integral is the probability that sitc-3 flipped 
up by a po-pi'ocess at /13 -I- 2| J| -I- h'" = 0. The second integrand is the probability that site-4 flips up at /14 + h" — 
and site-5 is down. Similarly the first integrand is the probability that site-2 fiips up at h' and site-1 is down. 
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Refer to figure (16) and the definition of P^i^r- Our object is to calculate the probability that site-3 flips down 
when site-2 and site-4 are up. This happens only if site-3 flips up after site-1 and site-5, i.e. sites 2, 3, 4 form a 
string of three down spins bordered by up spins at 1 and 5 just before 3 flips up. Because sites 2 and 4 are adjacent 
to up spins the distribution of /12 and /14 is the a posteriori (/'(/12) and 4>{hi) respectively. Also /13 > /i2 — 2\J\ and 
/13 > /14 — 2\J\ because 3 flips up before 2 and 4. We assume /14 > /12 and multiply our result by a factor of 2 to 
include the case /14 < /i2. Thus — 2\J\ < /13 < /12, /14 — /12 < 2\J\, and /12 < /14 < /12 + 2| J|. The probability that 
site-3 flips down when sitc-2 flips up is given by 

/ha rh' p-h'+2\J\ 

dh'4>{-h' + 2\J\) / dh"4){-h" + 2\J\) / d/i"'0(-/i"'-2|J|) (80) 

^00 Jh'-2\J\ J-h" 

The integrand of the first and the second integral is the probability that site-2 and sitc-4 flip up by a p2-proccss at 
h' and h" respectively. The last integrand is the probability that sitc-3 flips up by a pq- process at h'" . 

S. Ptits 

In figure(17) site-2 or sitc-4 could flip up first. Let us consider the case when site-2 flips up first i.e. h2> h^ — 2\J\. 
P-\itS is zero unless h2 < < /14 and /14 — 2| J| < h2 < h^. Therefore the probability that a new singlet is created at 
site-3 when site-4 flips up is given by 

pha ph' i.-h'+2\J\ 

Pnts{ha)^2 dh'4>i-h' + 2\J\) d/i"0(-/i")Pu(U |2;;/i") / dh"'<f>i-h"' -2\J\) (8I) 

J-00 Jh'~2\J\ J-h" 

The last integrand has the same interpretation as in the previous object P^i^r. The second integrand is the 
probability that site-1 is down when site-2 flips up by a pi-process at h". Finally the first integrand is the probability 
that site-4 flips up by a p2-pi'ocess at h' 

T. -Pf-J-tT 

Site-4 would flip up before site-2 in figure(17) if hi > h2 + 2\J\. Clearly h^> h2- Site-3 would flip down when site-2 
flips up if < h2 + 2\J\. Thus /13 and /14 lie in the range /14 - 2\J\ < < h2 + 2\J\ and /12 + 2\J\ < /14 < h2 + 4| J|. 
The probability that sitc-3 flips down when sitc-2 flips up is given by 

PnMha) = 2 / " dh'cj){-h')P^^{l i \2 i; h') ( dh"4>{-h" + 2\J\) [ ^ ' " dh"'cj){-h"' - 2\J\) (82) 

J-00 Jh'-2\J\ J-h" 

The integrand is similar to that in P-fits except that the first integral is for /i2 and the second for /14. 

U. Pf^^c/ 

P't'fW is associated with the destruction of objects characterized by P-^^-^q. The destruction can be analyzed on 
similar lines as their creation except that we now have the inequalities h2 < hi < h^ and ^13 — 2| J| < /12 < /i3. The 
probability that site-3 in P-\i-\q flips for the third time is given by 

Pmuiha) ^2 r dh'^{-h' + 2\J\) [ dh"<f>i-h")P^^{l i \2 h") f dh'" ^{^h"')P^^{5 i |4 U h'") (83) 

J-oo Jh'-2\J\ Jh'-2\J\ 

The second and the third integrands give the probability that sites 2 and 4 flip up at h" and h'" respectively. 
The first integrand is the probability that site-3 flips up by a p2-process at h' . Note that a priori distributions 
(/)(ft-2), 0(^3)7 0(^4) are used here because these sites have remained screened from site-1 and site-5. 
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P-\-\-^v is associated with the destruction of P^i^r. The inequahties that govern P-^^^y are h2 < < /?,3 + 2| J| and 
hs < h2 < hz -\- 2\J\. Thus the probabihty that the down spin in P-\i^r flips for the third time is 

dh'(t){~h' + 2\J\) / dh" ^{-h" + 2\J\) / dh'" 4){-h"' + 2\J\) (84) 

-oo Jh'-2\J\ Jh'-2\J\ 

The second and the third integrands give the probability that site-2 and site-4 flip up by a p2-pi'0cess at h" and 
h'" respectively. The first integrand is the probability that site-3 flips up by a p2-pi'ocess at h' . 



W. 

Pf^fw is associated with the destruction of -Pf4,ts- The appropriate inequalities for P-f^^w arc ^,3 < /14 < /13 + 2| J| 
and — 2| J| < /i2 < ^3- The probability that the down spin in question flips for the third time is 

Pmwiha) = 2 / ° dh'c^{-h' + 2\J\) [ dh"^{-h" + 2\J\) [ dh'"cf>{-h"')Pu{l i \2 I; h'") (85) 

J-oo Jh'-2\.]\ Jh'-2\.]\ 

The first integrand is the probability that site-3 flips up by a p2-pi'ocess at h' . The second and third integrands 
account for site-4 and site-2 respectively. 

X. P-ftt-f 

^ttt-f gives the fraction of singlets associated with P'^i^t that are destroyed at ha- These may be calculated 
in a similar manner as in the preceding case. We now have the inequalities /12 + 2| J| < /14 < /13 + 2| J| and 
/13 — 2| J| < /i2 < /13. Therefore the probability that the down spin in P-\i^^T flips up is given by, 

Pmxiha) = 2 / ° dh'(^{-h' + 2\J\) [ dh"<P{-h")P^^{l ; |2 h") f dh"'^{^h"' + 2\J\) (86) 

J-00 Jh'-2\J\ Jh'-2\.J\ 

The last two integrands pertain to sites 2 and 4 respectively, and the first to site-3. 



VII. MAGNETIZATION ON LOWER HYSTERESIS LOOP 



We are now in a position to write the magnetization on the lower half of the hysteresis loop. 



m{ha) = 1 - 2P^{ha) (87) 

where Pi{ha) is the probability that a randomly chosen site on the chain is down at applied field ha- A randomly 
chosen site on the chain can be characterized by the number n of up neighbors it has (n = 0, 1, 2). Thus we can write, 



Piiha) = PuM + Pmiha) + PuiiK) + Pmiha) (88) 

The first term corresponds to n 0, the next two terms that are equal by symmetry correspond to n = 1, and the 
last term corresponds to n = 2. We obtained the first three terms on the right-hand-side with relative ease in section 
(V). Surprisingly the evaluation of the last term i.e. P^^'[{ha) proved rather tedious requiring the calculation of 24 
terms as a pre-requisite [P-f^-|-yi(/ia) to P^^^x{ha)\- We have (so far) not found a simpler method to calculate P^x,-\{ha) 
in spite of much effort and thought. Putting all terms together we get. 
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PM = PmiK) + Pmiha) + Puiiha) + PmA{ha) + PnwiK) 



+-PntQ('^a) + P^i^RiK) + PntsiK) + PntTi^a) 



(89) 



We note that several objects that we calculated in the preceding section do not appear explicitly in the above 
equation e.g. P-xi^siha) does not appear explicitly in equation (89). However it is necessary to calculate P-fi-fsiha) 
because it is needed in the calculation of P-fi-\-M{ha) and P-fi-^oiha) that appear in the final formula. Similar remarks 
apply to other terms that were calculated but do not appear explicitly in equation (89). 

The magnetization in increasing field is obtained by substituting equation (89) in equation (87). The magnetization 
in decreasing field on the upper half of the hysteresis loop is given by symmetry, 



In the next section we compare the theoretical result against numerical simulations of the model in selected cases. 
We consider a uniform bounded distribution of the quenched field as well as a Gaussian distribution. As may be 
anticipated, the agreement between theory and numerical simulation is quite good. 



Simulations of the model have played an important role in the analysis presented here. Although an exact analytic 
result has to be necessarily in agreement with the simulations within numerical errors but arguments based on 
conditional probabilities can be subtle and prone to errors. Therefore at each step of the analysis, we devised a 
simulation of the model to yield the probability of the event being calculated. Occasionally the two would not match 
in the first instance necessitating a rethink of the analysis and locating the error in the argument. Thus each of the 
theoretical expression in the preceding sections was verified by simulation of the model for a bounded distribution of 
the random field with half-width A = 1.25|J| and a Gaussian distribution with standard deviation a = .5\J\. The 
comparison between the theoretical hysteresis loops and those obtained by simulation was shown already in section 
III. Here we show the comparison for a few other quantities that enter the expression for the hysteresis loops. 

We begin with the probability per site of a doublet on the lower half of the hysteresis loop. Figure (18) shows 
the theoretical expression for P4.4. (/la) for the Gaussian distribution with the corresponding data from numerical 
simulation superimposed on it. The fit is so close that the two are indistinguishable on the scale of the figure. Figure 
(19) shows similar comparison for a uniform distribution with A ~ 1.25. In each of the following figures data from 
the corresponding numerical simulation has been plotted along side the theoretical expression. Even numbered figures 
are for the Gaussian distribution and the odd numbered figures for the uniform distribution. In some cases the match 
between theory and simulation is so good that there appears to be only a single curve in the figure. In other cases 
( when the probability of the event is relatively small and finite size corrections are larger) we can barely make out 
that there are two curves that almost lie on each other. Figure(20) and figure (21) show the result for Pin,{ha) for 
the Gaussian and the uniform distribution respectively. Results for P-\u{ha) + Pu^iha) are shown in figures (22) 
and (23) for Gaussian and uniform distributions respectively. Figure (24) shows the comparison between theory and 
simulation for P-^^-^Aiha) in the case of Gaussian distribution. Figure (25) shows similar comparison for P^i^oiha) 
for the uniform distribution. Figure (26) is for Pf^fAiiha) for the Gaussian distribution. Figure (27) is for P-ft-fAf (^a) 
for the rectangular distribution. Figures (28) and (29) each contain two objects. Figure (28) shows P^i^giha) and 
P-\^^uiha) for Gaussian distribution. Figure (29) shows P-\i^B,{ha) and P-\-^-\-v{ha) for uniform distribution. 

In conclusion we have obtained the zero-temperature hysteresis loop of a one dimensional anti-ferromagnetic random 
field Ising model in the case when the driving field varies from —00 to 00 and back to —00 infinitely slowly. The 
problem is simple to state but difficult to solve. The theoretical result for the hysteresis loop involves integrals that 
have to be evaluated numerically in most cases. We have shown that our results fit numerical simulations of the 
model quite well. However, several aspects of the problem still remain unsolved. For example we have obtained the 
hysteresis loop when the driving field takes the system from one saturated state {ha = —00) to another {ha = 00). 
In this case we have a complete knowledge of the statistical history of the system i.e. if a site is up at ha, we know 
the relative probability of different sequence of events that result in this site being up. We are not in a position (so 
far) to obtain the hysteretic response of the system starting from an arbitrary initial state. It would be interesting 
to have an analytic solution of the problem in higher dimensions as well. Numerical simulations suggest that the 



m"{ha) 



m{-ha) 



(90) 



VIII. COMPARISON WITH SIMULATIONS AND CONCLUDING REMARKS 



anti-ferromagnetic hysteresis loops in higher dimensions have several plateaus for low values of A and a as compared 
with I J|. Exact analytic solutions of problems with quenched disorder are uncommon in statistical mechanics. One 
may even ask if they are worth the effort that has to be put in trying to obtain them. However exact solutions are 
intellectually satisfying and provide a framework for understanding a wide class of complex phenomena. We hope 
there will be more progress in this direction in the future. 
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FIG. 1: Hysteresis loop for an anti-ferromagnetic random-field Ising model with J = —1 and A = 0.5 (see text). The a;-axis 
shows the applied field and the y-axis magnetization per spin. As |J| < A, each half of the hysteresis loop comprises three 
ramps separated by two plateaus. The lower half of the loop shows magnetization in increasing field and the upper half in 
decreasing field. A theoretical expression has been superimposed on the numerical data. 
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FIG. 2: A magnified view of the theoretical and simulation hysteresis loops for A — 0.5 where the t/-axis shows the magnetization 
in increasing and decreasing field as measured from the average of the magnetization on the lower and the upper half of the 
hysteresis loop in figure [T]at corresponding applied field. 
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FIG. 3: Hysteresis loop for an anti-ferromagnetic random-field Ising model with J = —1 and A = 1.25 (see text). The x-axis 
shows the applied field and the y-axis magnetization per spin. As |J| > A, the plateaus of figure [1] disappear and the three 
ramps merge into each other. A theoretical expression has been superimposed on the numerical data. 
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FIG. 4: Theoretical and simulation hysteresis loops for A — 1.25 where the magnetization along increasing and decreasing 
field is measured from the average of the magnetization on the lower and the upper half of the hysteresis loop in figure [3] at 
corresponding applied field. 
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FIG. 5: Hysteresis loop for an anti-ferromagnetic random-field Ising model with J = — 1 and a — 0.5 (see text). Notice the 
approximate similarity with the hysteresis loop in figure [1] but the absence of sharp ramps and plateaus. For a Gaussian 
distribution the three ramps merge into each other for any value of a although this is less pronounced at smaller values of a. 
A theoretical expression has been superimposed on the numerical data. 
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FIG. 6: Magnified theoretical and simulation hysteresis loops for a — 0.5 where the magnetization along increasing and 
decreasing field is measured from the average of the magnetization on the lower and the upper half of the hysteresis loop in 
figure [S]at corresponding applied field. 




FIG. 7: Spins on Ramp-I in an applied field — 2| J| — h. Filled circles show sites with quenched field hi > h. The probability per 
site of a doublet (two adjacent down spins) such as AB is equal to e^^*", where p is the fraction of filled circles on the infinite 
lattice. 



\ h2 




FIG. 8: A doublet on Plateau-I: hi and /i2 are the quenched random fields on the doublet sites 1 and 2 respectively. 
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FIG. 9: Two adjacent doublets on Plateau-I: Each doublet separates the lattice into two parts whose evolution histories on 
Ramp-I are independent of each other. Evolutions inside each dashed box is shielded from outside. The probability that spin 
at site 3 flips up on Ramp-I is therefore equal to ^ . Given this, the probability that the spins at sites 1 and 5 remain down all 
along Ramp-I is equal to i each. The shielding property of the boxes can also be used to determine a posteriori distribution 
of random fields hi, hi, h^, hi, and /15. 
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FIG. 10: A singlet (site 3) with one next nearest neighbor down (site 1), and one next nearest neighbor up (site 5). When the 
singlet turns up at an applied field ha, the spin at site 2 stays up if A < | J|, but the spin at site 4 fiips down if hi < hz. 
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FIG. 11: Two adjacent singlets on Plateau-I: If /12 ~ mm(/i2, ^4), and /13 < /12, then the spin at site 3 will flip down when the 
spin at site 2 flips up on ramp-III. This process creates a new singlet on ramp-III. 




FIG. 12: A singlet followed by a doublet on Plateau-I: If > h^, and /13 < /i2, then a new singlet will be created at site 3 
when the spin at site 2 turns up on ramp-III. 
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FIG. 13: A doublet followed by a singlet on Plateau-I: 
when the spin at site 5 turns up on ramp-III. 
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[26] 
[27] 
[28] 
[29] 
[30] 
[31] 
[32] 
[33] 
[34] 
[35] 



[36] 
[37] 
[38] 
[39] 
[40] 



X Ilia, M L Rosinberg, P Shukla, and E Vives, Phys Rev B 74, 224404 (2006). 

X lUa, M L Rosinberg, and G Tarjus, Eur Phys J B54, 355 (2006). 

P Ayyub et al, Phys Rev B 57, R5559 (1998). 

K Fukuma and M Torn, Earth Planets Space 50, 9 (1998). 

I Chiorescu et al, Phys Rev Lett 84, 3454 (2000). 

E E FuUerton et al, Appld Phys Lett 77, 3806(2000). 

O Waldmann et al, Phys Rev Lett 89, 246401(2002). 

K Takanashi, Appld Phys Lett 63, 1585 (1993). 

Khian-Hooi Chew et al, Appld Phys Lett 77, 2755 (2000). 

J Kisker,H Rieger, and H Schreckenberg, J Phys A: Math. Gen. 27, L853 (1994). This paper discusses the non-equilibrium 
dynamics of a non-random, one-dimensional Ising model, with three-spin interactions, at low temperatures. It is qualita- 
tively similar to the dynamics of the random field Ising model at zero-temperature. 
H Reiger, Physica A 224, 267 (1996). 

F Ritort and P SoUich, Advances in Physics 52:4, 219 (2003). 

C Toninelh, PhD thesis (2003). 

J W Evans, Rev Mod Phys 65, 1281 (1993). 

L G Mityushin, Prob Peredachi Inf 9, 81 (1973). Also see reference (39) for a detailed discussion of the screening property 
of this class of problems. 




FIG. 14: An up spin followed by a doublet. Site-1 may have flipped up under a po-process or a pi-process. Subsequently when 
site-3 flips up we get a singlet at site-2. The fraction of such singlets depends on the details of how sites 1 and 3 have flipped. 
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FIG. 15: A doublet followed by a doublet. We get a singlet at site-3 if sites 2 and 4 flip up before sites 1 and 5 and site-3 flips 
down because it is unstable with both neighbors up. 
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FIG. 16: A singlet followed by a singlet. A singlet is created at site-3 if sites 2 and 4 flip up under a p2-process and then site-3 
flips down because it is unstable when both neighbors are up. 




FIG. 18: Probability of a doublet P^^(/ia) for a Gaussian distribution with a — 0.5| J|. Simulation data has been superimposed 
on the theoretical expression. 
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FIG. 20: Probability of Piii{ha) for a Gaussian distribution with a — 0.5|J|. Simulation data has been superimposed on the 
theoretical expression. 
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FIG. 21: Probability of Piii{ha) for a uniform distribution with A — 
theoretical expression. 
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FIG. 22: P-fii{ha) + Pii^{ha) for a Gaussian distribution with a = 0.5|J|. Simulation data has been superimposed on the 
theoretical expression. 
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FIG. 23: P^ii{ha) + Pii^{ha) for a uniform distribution with A = 1.25| J|. Simulation data has been superimposed on the 
theoretical expression. 



□ .3 



Q.2 ■ 



□ .1 ■ 



□ 




-3 -2 -1 



FIG. 24: Theory and simulation for P^i^Aiha) for a Gaussian distribution with a = 0.5| J|. 
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FIG. 26: Theory and simulation for P^^^Ai{ha) for a Gaussian distribution with a = 0.5|J|. 
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FIG. 27: Theory and simulation for PfffN{ha) for a uniform distribution with A = 1.25|J|. 
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FIG. 28: Theory and simulation for P-iitQi^a) on the left, and P-iffu{ha) on the right for a Gaussian distribution with 
a = 0.5|J|. 



